Incoherent Transport through Molecules on Silicon in the vicinity of a Dangling Bond 
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We theoretically study the effect of a localized unpaired dangling bond (DB) on occupied molec- 
ular orbital conduction through a styrene molecule bonded to a n^"*" H:Si(001)-(2xl) surface. For 
molecules relatively far from the DB, we find good agreement with the reported experiment using 
a model that accounts for the electrostatic contribution of the DB, provided we include some de- 
phasing due to low lying phonon modes. However, for molecules within 10 A to the DB, we have to 
include electronic contribution as well along with higher dephasing to explain the transport features. 
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I. INTRODUCTION 

Interest in atomic scale conduction continues to grow 
at a rapid pace both in the engineering and scientific com- 
munities. The comparatively stable C-Si bond and the 
well characterized Si surface chemistrji^ not only place 
organic molecules in an ideal position to complement ex- 
isting Si technology, but also offer a unique platform to 
study scientific phcnomcn^i^ with experimental repro- 
ducibility. Notably, Si contacts may give rise to novel 
features such as negative differential resistance^i^i^, al- 
though the mechanism is not yet established^. 

In Si technology, dangling bond (DB) defects in the 
form of Pb centers have created reliability issues - con- 
sider, negative bias temperature instability in metal- 
oxide-Si structures. However, one could engineer such 
defects to be useful in novel devices. Recently Piva et 
al^ observed significant shifts in the onset voltages along 
a styrene chain as function of distance from a negatively 
charged unpaired DB on a n++ H:Si(001)-(2x 1) surface 
using scanning tunneling spectroscopy (STS). In this pa- 
per, we theoretically study their results. A comparison 
between calculated and experimentally observed current- 
voltage (I-V) plots is shown in Fig. 1. The DB, being 
negatively charged, pushes molecular levels upward. A 
schematic energy level/band diagram is shown in Fig. 
1. The onset voltage depends on the gap between the 
equilibrium chemical potential (fio) and molecular levels. 
Given that the onset voltage is decreasing in the vicinity 
of the DB, one concludes that conduction is through oc- 
cupied orbitals. The trend would have been reversed for 
unoccupied orbital conduction. 

Furthermore, observed onset voltages for molecules 
AA and 8A^ away from the DB are grouped together; 
whereas a molecule 12A away has a distinctly different 
onset voltage. Electrostatically, one would expect them 
to follow an inverse distance relationship. To explain the 
above discrepancy, we propose that the DB also affects 



the I-V characteristics of the styrene chain electronically 
by introducing a near-midgap state in the local density 
of states (LDOS) of Si atoms AA and 8A away. For these 
molecules, this contribution would lead to an early onset 
of conduction at approximately 0.5V (rather than I.IV 
due to the bulk Si band gap) as shown in Fig. 1. Since 
the DB wavefunction decays exponentially, this contri- 
bution would be nearly absent for a molecule 12 A away. 
Apart from this, the effect on transport due to the DB 
state would be pronounced when the highest occupied 
molecular orbital (HOMO) level is sufficiently broadened 
to form a significant overlap with the DB state as shown 
in the inset of Fig. 2. Moreover, this effect would be 
diminished if the HOMO level is too far from the DB 
state. 



The impact of dephasing due to low lying phonons on 
experimentally observed lineshapes is further explored in 
this paper. In molecular junctions, the dominant phonon 
modes are the ones having low energy, which arise due 
to the combined motion [translational, rotational, etc] 
of the entire molecule relative to the contacts and have 
been discussed theoretically in RefsJ^ii^. Experimen- 
tally, some evidence of the presence of these low lying 
modes is availablei^iiiii^. Since experimental results are 
suggestive, a detailed study needs to be done to establish 
it on a solid footing. These low lying phonon modes are 
excited due to inelastic electron-phonon scattering. Sim- 
ilarly, the styrene should be able to vibrate, rotate and 
move with respect to the Si and tip contacts, and hence 
we expect a similar trend. Since hw <^ ksT at room tem- 
perature for these modes, we include their effect through 
elastic dephasing for simplicity. Furthermore, broaden- 
ing due to such modes could be large as in Ref.— . The 
dephasing strength used in our calculation gives a com- 
parable broadening. 
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FIG. 1: (color online) Calculated transport properties of a 
styrene molecule bonded to a n^""" H:Si(001)-(2x 1) surface in 
the vicinity of a negatively charged DB. Experimental data is 
indicated by x marks^ii^. As a function of decreasing distance 
from the DB, the onset voltages shift toward lower tip voltages 
implying occupied molecular orbital conduction. I-V plots 
for molecules and 8A away have grouped onset voltages. 
The electrostatic and electronic contributions of the DB are 
included in these calculations. We assume a lower dephasing 
strength (D) of lAeV^ for the molecule that is farthest from 
the DB, as compared to 1.85eV^ for the closer ones. 



II. THEORETICAL MODEL AND 
ASSUMPTIONS 

We use the single particle non-equilibrium Green's 
function formalism in the mean field approximation using 
a non-orthogonal basis to model quantum transport as in 
Refjii. We extend this work by including elastic dephas- 
ing within the self-consistent Born approximatio n^^i^^ . 
The time-retarted Green's function is defined as: 



G=[iE + iO+)S -Hd-Ud-UoB-^ 



Si, tip 



where is the device Hamiltonian and S is the overlap 
matrix. Ud is the potential experienced by the molecule 
and consists of (1) Laplace potential due to applied tip 
voltage, (2) Band lineup potential due to the Fermi level 
mismatch of n++-Si and tip contacts, (3) Hartree and 
Image potentials, due to non-equilibrium statistics of the 
electrons on the molecular region. Udb is the poten- 
tial due to the negative charge on the unpaired dangling 
bond, which causes a shift in the molecular levels. The 



electrostatic and electronic boundary conditions are set 
by the tip voltage and contact self-energies {^si,tip) 
spectively. Ss is the scattering self-energy. In a non- 
orthogonal basis for the elastic dephasing, the scattering 
broadening function is given as2£: 



D SA{E)S 



(2) 



where S™ is the scattering inflow function and 1]°"* is the 
scattering outflow function. D is the dephasing strength 
which is a fourth rank tensor, but in our calculations it 
is approximated by a scalar and is used as a fitting pa- 
rameter. Treating D as a scalar would ignore the minute 
details of dephasing that may be otherwise present when 
a carrier gets scattered from one electronic state to an- 
other electronic state more strongly or weakly. A is the 
Spectral function defined as i{G — G^) = G" + Gp, where 
G" and G^ are the electron and hole correlation functions 
respectively. Using Eq. 2, the scattering self-energy is 
given as: 
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It has been shown^i that the real part of ^^(i?) [cal- 
culated by the Hilbert transform as shown in Eq. 2] is 
about few meV for different organic molecule contacted 
by Au(lll) electrodes due to slowly varying TsiE). We 
expect it to be small for the molecular transport through 
Si valence band due to slowly varying density of valence 
band states. However, this approximation may not be 
appropriate for a one dimensional contact like carbon 
nanotube or graphene nanoribbon due to van Hove sin- 
gularities. Given that this calculation is computationally 
very expensive and has small effect on the final results, we 
ignore the real part of T.,{E). Since S^, G, G"(== ~iG<) 
and Ss"(= — «E<) depend on each other, we solve for 
the above four quantities self-consistently along with the 
Hatree self-consistent loop of the applied tip voltage. 

The electronic structure of Si(OOl)^, styrene^ and the 
tip is calculated using Extended Hiickel Theory (EHT). 
Importantly, EHT provides the correct bulk Si band gap 
(as compared to LDAi) and captures the surface and 
band structure effects accurately^^i^. We incorporate 
an isolated DB in a H:Si(001)-(2x 1) symmetric dimer 
unit cell, composed of 64 Si atoms with 8 at the surface, 
out of which 7 are passivated and one is unpassivated 
which has the dangling bond, as shown in the inset of 
Fig. 4 (visualized using Gauss View^l) . Additional sur- 
face relaxation due to the DB is small-^ and hence ig- 
nored. We calculate the surface Green's function (gs) 
recursively^^ and find that the DB gives rise to a near- 
midgap state in the Si band gap similar to Refi^. We 
assume a constant gs for the tip similar to Refi^. How- 
ever, Etip is still energy dependent, since it is defined as 
[{E + iO+)Si - Hi]gs[{E + iO+)Sl - where Si and 
Hi are the overlap and Hamiltonian matrices between tip 
and styrene molecule. For a good STS tip which can give 
atomic resolution, the last tip atom at the apex domi- 
nates STS and hence we use a single tungsten (W) atom 
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as shown in Fig. 2(a). This apparent barrier height is 
about lOeV with original parameters. Apart from this, 
the HOMO level of W atom hes at -10.73eV, whereas 
the work function of W(llO) tip is around 5.2eV, there- 
fore we have an offset of 5.53eV with respect to vac- 
uum. Similarly, the Si conduction band edge (Ec) ob- 
tained from EHT has an offset of 7.75cV. To calculate 
the offset of styrcne, we use the HOMO level obtained 
from DFT-B3PW91/6-311g* calculations^ as a refer- 
ence, since the occupied states in DFT are relatively well 
characterized^. To ensure a proper alignment with re- 
spect to vacuum, we shift the Hamiltonians by the above- 
mentioned offsets Eof fset-i,j for the systems i and j 
within EHT procedure as follows: 
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FIG. 2: (color online) Tip modeling and electrostatic contri- 
bution of the DB. (a) Tunnel current as a function of tip height 
showing exponential decay, which gives an apparent barrier 
height of 4.3eV. (b) A 3D continuum finite-element electro- 
static calculation showing a hyperbolic trend in the average 
potential, due to a negatively charged DB, along a styrene 
chain. The insets show a contour plot potential profile 4A 
away. 



as tip. Transport quantities are then calculated through 
a single molecule. The molecular structure is taken from 
Refi^, which reports structure for a styrene molecule on 
hydrogenated Si surface. There is an important difference 
between bonding geometry of styrene on hydrogenated 
and unpassivated Si. For hydrogenated Si, one C atom 
of styrene bonds to one Si atom as shown in Fig. 1. For 
unpassivated Si, styrene goes through a cyclo-addition 
process and the two C atoms on the styrene bind with 
two atoms of a Si dimer [not shown in this paper] . 

For Ylfip, the EHT parameters of the s-orbital basis set 
arc modified^S (similar to Refi^) to get an appropriate 
apparent barrier height (0ap)- We obtain (pap = 4:.3eV 
from the calculated l(z) plot^ with 2V applied at tip 



^ij ~ 2^ij{KiHii+KjHjj+Eoffset-i+Eoffset~i) (4) 

In order to match the experimentally observed onset 
voltages, we shift the molecular levels by -|~0.5cV [shifted 
toward the valence band minimum]. This shift is a one 
time adjustment for all the molecules under considera- 
tion and hence it affects the onset voltages of all the 
molecules by an offset of -F0.5V. Further changes in the 
onset voltages come due to the DB potential Ud b , which 
varies inversely with distance of different molecules from 
the DB. Since, the density of states inside HOMO-LUMO 
(lowest unoccupied molecular orbital) gap is quite small, 
a shift of 0.5V can be caused by a very small change 
in Hartree potential at equilibrium due to small charge 
transfer. However, the lineshapes of the calculated I-V 
characteristics are not sensitive to this shift and they re- 
main convex [see Sec. HI]. Incorporating this shift results 
in the HOMO level being about 0.4eV below the valence 
band edge [after incorporating the effect of Si contact as 
well in the form of Si sclf-encrgy - the real part of the 
self-energy shifts the molecular levels]. This means that 
HOMO level is about 0.9eV below the valence band edge 
before applying this shift. The position of HOMO level 
in our study is different from Rcf^^, where it is reported 
to be about IcV below the valence band edge. How- 
ever, the difference in our case is that the styrene chain 
is in the vicinity of a charged dangling bond which could 
change local nature of Si-C bond. This governs the charge 
transfer at equilibrium, whereas in Refi^, there is no dan- 
gling bond. Furthermore, we do not calculate the charge 
transfer at equilibrium, rather use the shift in energy lev- 
els due to charge transfer as an adjusting parameter. In 
short; there are two fitting parameters in these calcula- 
tions, shift of energy levels [a one time adjustment] and 
dephasing strength D. These two fitting parameters are 
used to capture two independent features in the experi- 
ment. The first one, i.e. shift of energy levels, is used to 
quantitatively match the experimentally observed onset 
voltages and it does not affect the lineshapes observed 
whether we include dephasing or not. The second pa- 
rameter, i.e. dephasing strength D, is used to reproduce 
the experimentally observed concave lineshapes. The de- 
phasing affects the onset voltage in a minor fashion since 
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FIG. 3: (color online) Effect of dephasing and electrostatic 
contribution of the DB on transport properties, (a) Calcu- 
lated I-V plots for the molecule 12 A away with and with- 
out dephasing. Without dephasing, the lineshape [red (dark 
gray) dotted line] does not match with the experiment (blue 
X marks). With D = lAeV^ , there is a good quantitative 
agreement with the experiment, (b) Calculated I-V plots for 
molecules AA, SA and 12A away from the DB after incorpo- 
rating Udb- 



it broadens the energy levels and hence current starts 
flowing at lower tip voltages. 

We solve the 3D Laplace's equation using the finite- 
element method to obtain the potential profile across 
styrene due to the tip voltage. This method has been 
discussed in detail in Ref.— . This method also provides 
the zero bias band lineup potential due to the Fermi 
level mismatch of l.lSeV between the tip and r7,++ Si. 
Solving for a 3D potential ensures that Stark effect is 
captured for the applied tip voltage. Furthermore, for 
the experimental conditions the molecular levels shift by 
about one tenth of the applied tip voltage and/or Fermi 
level mismatch, since the tip is quite far from the styrene 
molecules. This scenario is different from the ones re- 
ported in Refsi^i^, where the molecular levels shift by 



about one third of the tip voltage because tip is quite 
close to the molecule. Apart from this, the tip heights 
are modified^ to reflect actual changes during experi- 
mental data acquisition^ and hence tip heights are not 
used as adjustment parameters. Rather, it implies that 
our model is flexible to incorporate the very delicate ex- 
perimental details. 

For calculating the potential due to the DB, a 3D finite- 
element continuum calculation^ is performed by taking 
into account the DB's approximate shapoSI. The styrene 
chain is approximated by a lOOA long 3D box with width 
and height corresponding to that of styrene and having 
a relative dielectric constant of about 2.4. This box is 
placed on a Si bulk having a relative dielectric constant 
of about 12. An STM tip is placed about Inm above 
the styrene chain in the form of a metallic sheet. An 
average of the calculated potential {Udb) taken over the 
lateral cross-section is included in our transport calcula- 
tion. A plot of Udb as a function of distance from the 
DB is shown in Fig. 2(b). The potential profile across 
a molecule away from the DB [shown in the inset of 
Fig. 2(b)], shows a weak electrostatic contribution in Si. 
Incorporating Udb in an averaged fashion ignores Stark 
effect within a single molecule due to charge on the dan- 
gling bond. This may quantitatively affect the results in 
a minor way because as seen in the inset of Fig. 2(b), the 
dangling bond potential is comparatively constant over a 
single molecule, i.e within about 0.5-0.7eV. This simple 
continuum model seems to be satisfactory for calculating 
the electrostatics of the DB, because (1) the screening 
effect of the styrene chain should be small due to the one 
dimensional nature of the styrene chain. In particular, 
occupied states are more localized than unoccupied states 
and hence this screening is expected to be small for occu- 
pied states conduction [a case discussed in this paper] (2) 
we obtain a hyperbolic dependence as a function of dis- 
tance from the dangling bond i.e. Udb ~ 2.4(ey)/d(A), 
which is expected. For calculating the electrostatic con- 
tribution, treating the charged DB, styrene chain and 
Ti^"*" doped Si in atomistic manner may result in quan- 
titatively different Udb, however it is still expected to 
show a hyperbolic dependence and can be readily incor- 
porated in our model. Such a change would result in 
different onset voltages as a function of distance from 
the dangling bond, however the reported lineshapes in 
Sec. Ill would remain the same. 

Moreover, based on a finite-element MEDICS^ cal- 
culation, we conclude that tip induced band bending 
is also negligible. Additionally, dopants in n~^~^ Si in- 
troduce states approximately 50meV— below Ec, which 
contribute weakly^ to transport for a tip voltage of 
about 50meV. In our calculations, we ignore the above- 
mentioned effects. The Hartrce potential for the molecule 
is included via the complete neglect of differential over- 
lap method^. Image effects are incorporatedii to ensure 
that self-consistency is not over-estimated, which can al- 
ter results significantly. In all of the above electrostatic 
calculations, the tip is taken to be a metal sheet having 
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FIG. 4: (color online) Electronic contribution of the DB. The 
calculated local density of states of Si atoms, 4A and 8A from 
the DB, shows that the DB introduces a near-midgap state. 
Inset shows the surface portion of the unit cell used. 
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III. DISCUSSION OF RESULTS 



A transport calculation to study the effect of dephasing 
for a molecule 12A away from the DB is presented in Fig. 
3(a). Without dephasing, the lineshape [red (dark gray) 
dotted line] does not match with the experiment [blue 
(dark gray) x marks]. The lineshapes for the theoret- 
ically calculated and experimentally observed I-V char- 
acteristics are quite different. Experimentally observed 
lineshapes are concave, whereas the calculated lineshape 
has a convex character. For coherent transport, obtain- 
ing a convex lineshape is expected because it represents 
that current starts flowing through a level and then satu- 
rates when the level is fully inside the bias window. The 
corresponding conductance peak [not shown] has a full 
width half maximum (FWHM) of about 0.2eV. This is 
roughly equal to the coupling between styrene and Si con- 
tact because tip is quite far and hence the TV character- 
istics do not get affected much by the Hartree potential 
due to the non-equilibrium carrier statistics. This small 
coupling is expected since styrene binds with Si through 
only one Si-C bond and HOMO is localized on the aro- 
matic ring of the styrene molecule, see e.g. Refs.— i^. 
However, after including dephasing (with D = l.AeV'^), 
not only is the lineshape reproduced, but good quantita- 
tive agreement with the experiment is achieved. 

This transition from a convex lineshape to concave 
lineshape needs further discussion. The physical effect 
of dephasing is that it broadens the molecular levels. 



FIG. 5: (color online) Effect of electronic and electro- 
static contributions of the DB on transport properties, (a) 
T{E,V)[ftip — fsi] plots (where T(E,V) is the transmission 
and flip, Si are the Fermi functions for the contacts) for 
molecules 4A, 8A and 12 A away from the DB. For molecules 
4A and 8A away, there is conduction inside band gap due 
to electronic contribution of the DB. (b) Corresponding I-V 
plots. Onset voltages for molecules 4 and 8A away match well 
with the experiment, whereas lineshapes do not. 



Since the local density of states (LDOS) of Si contact 
[see Fig. 4] increases with decreasing energy /increasing 
positive tip voltage, the current starts increasing in a 
fashion proportional to the coupling to the Si contact, 
which is proportional to the LDOS of Si contact. In 
this case, as can be seen in Fig. 4, the LDOS increases 
with decreasing energy as one goes away from the va- 
lence band edge. The current increases with increasing 
tip voltage in a concave manner because of this increase 
in LDOS of Si contact. This LDOS corresponds to in- 
creased surface Green's function of Si and hence larger 
Si self-energy Sgi, which results in a higher current. In 
order to adjust for this, we have to use higher tip height 
of 8.67 A in our model, after including dephasing, instead 
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of 7.79A for coherent transport to get the same current 
level. Since experimentally tip heights in general remain 
to be unknown, this variation is acceptable. However, 
relative tip heights are known experimentally. In the set 
of experiments we consider, these relative heights were 
changed^ during experimental data acquisition for the 
molecules 4A and 8A away from the DB and wc include 
them in our model24. 

Using the same D for molecules 4A and SA away from 
the DB, and including Udb, we obtain the I-V plots 
shown in Fig. 3(b). Although the lincshapcs are similar 
to the experimentally observed lineshapes, there is still 
a significant mismatch in the onset voltages between ex- 
periment and calculation for molecules 4A and 8A away. 
Experimentally, their onset voltages are grouped together 
and are much smaller than the calculated ones. This is 
expected since the valence band edge corresponds to a 
tip voltage of about I.IV and hence very small current 
should be expected for tip voltages below I.IV. However, 
as can be seen in Fig. 3(b), experimentally conduction 
starts from about 0.5V. We interpret this disagreement 
between experiment and calculation to be induced by the 
DB. Besides contributing electrostatically, the DB intro- 
duces a near-midgap state in the LDOS of Si atoms be- 
neath the styrcne chain extending up to be about lOA. 
The calculated LDOS for Si atoms AA and 8A away from 
the DB is shown in Fig. 4. A near-midgap state is clearly 
visible inside the band gap which is responsible for the 
early onset of conduction at about 0.5V of tip voltage. 
The local DOS contribution of the DB state exponen- 
tially decays and is nearly absent for a styrene molecule 
about I2A away. 

In our transport model, the DB is taken to be fully oc- 
cupied under equilibrium and non-equilibrium conditions 
through inelastic processes inside the Si contact. This is 
a valid assumption since the energetic location of the DB 
state is well below the chemical potential. Furthermore, 
the electronic structure calculations are for intrinsic Si, 
where the effect of dopants and charged DB is ignored 
for simplicity, otherwise the calculations become com- 
putationally prohibitive. The effect of doping, however, 
is electrostatically included in our transport calculations 
by adjusting the equilibrium Si chemical potential. Be- 
sides this, treating the DB state as charged is expected 
to cause a shift in its energetic position within band gap. 
Since these states arc thermally broadened on the order 
of 0.5cV or mora^ at room temperature and that the 
DB state is close to midgap, the results presented are 
expected to be comparatively insensitive to this shift. 

In Fig. 5(a), the bias dependent calculated transmis- 
sion T(E,V) through molecules AA, 8A and 12A away 
from the DB is shown. For molecules 4A and 8A away, 
there is finite conduction inside the band gap due to this 
DB state, which is absent for the molecule 12A away. 
For the molecule iA away, at higher tip voltages, the 
level below the HOMO level i.e. HOMO-1 level, also 
starts conducting. Fig. 5(b) shows the corresponding I- 
V characteristics after including the DB's electronic and 



electrostatic contribution with D = l.ieV'^. The onset 
voltages for molecules AA and 8 A away match reasonably 
well with the experiment and are grouped together, but 
the lineshapes still differ. There is an increase in current 
at about 0.6V and then current tends to saturate. This 
is because contribution of the DB state is small in this 
energy range. Using a higher D of 1.85eV^ for molecules 
AA and 8A away broadens the molecular spectral density 
more and hence reproduces the experimentally observed 
lineshapes as shown in Fig. 1. These molecules, being 
at the end of styrene chain, may have more phonon de- 
grees of freedom thus resulting in a higher dephasing. 
Another possible explanation could be that the dephas- 
ing processes inside the Si contact broaden the LDOS 
contribution of the midgap state. Further experimental 
and theoretical work needs to be carried out to estab- 
lish a detailed understanding of the dephasing mecha- 
nism. For example, conducting temperature dependent 
transport experiments would sharpen the spectroscopic 
features with decreasing temperature, since D oc T for 
hoj <C fc^T. Apart form this, using different molecules 
with varying HOMO levels may give an insight into the 
electronic contribution. 

Since we are calculating transport properties of an iso- 
lated styrene molecule and comparing with the styrene 
molecules embedded in a chain, the question arises if we 
arc missing something. As far as one can obtain atomic 
resolution with the tip which means that there is one 
apex atom and the tip is above a styrene molecule, the 
conduction from the tip to other styrene molecules in- 
side the chain is expected to be small. Because styrene 
molecules are about AA away and if the tip is about 
8A above the styrene molecule, the distance of the tip 
atom from the next styrene molecule would be about 9 A. 
This should result in an order of magnitude lower current 
through the next styrene molecule. Therefore, although 
the conduction along the styrene chain may be notice- 
able as discussed in Ref4, since the tip is placed above 
a styrene molecule, the lowest conduction path is still 
through this particular styrene molecule and hence con- 
duction through the neighboring styrcne molecules would 
still be small. However, if the tip is placed between two 
styrene molecule [a case not discussed in this paper], 
one has to consider transport though both the styrene 
molecules as well as tunneling between them. Since in 
our case, the neighboring styrene molecules are expected 
to remain in equilibrium, the effect of self-consistent po- 
tential variation in the neighboring styrene molecules can 
be ignored as well. It should be noted that our case is 
different from that of a self-assembled monolayer device 
contacted by a flat contact, because in this case all the 
molecules are expected to conduct and hence a molecule 
would be affected by the self-consistent potential of the 
neighboring molecules. The self-consistent potential of 
the styrene molecule through which the conduction is 
occurring could affect the neighboring styrene molecule 
electrostatics and in return the perturbed neighboring 
molecules could affect the self-consistent potential of the 
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conducting molecule itself. Since, tip is far from the 
molecule and self-consistent potential is small [although 
not negligible], this effect is expected to be small. 

We also ignore the tilt angle of the styrene molecule 
with respect to the Si substrate as a function of distance 
from the dangling bond that may be present in reality. 
It has been shown previously that the tilt angle for a 
styrene chain in the absence of the dangling bond as a 
function of distance from the end of the styrene chain 
may lead to significant changes in the current flowing 
through the styrene chain and hence the observed STM 
images. This end effect has been observed experimentally 
for low bias and qualitatively reproduced theoretically 
in Refj^. This end effect vanishes at higher bias and is 
still not well understood. As shown experimentally in 
Refi^, in the presence of a charged dangling bond at the 
end of the styrene chain, this end effect is not present. 
Therefore, assuming that the tilt angle of the molecules 
near the end of the chain does not change seem to be 
a good approximation. However, a detailed quantitative 
study is needed. It would also be interesting to analyze 
how the low lying phonon spectrum gets modified for a 
styrene chain on Si surface in the presence of a dangling 
bond, whereas previously phonon spectrum of a single 
molecule on a substrate has been discusse d^^'^^ . The 
inter-molecular interactions are expected to give rise to 
new peaks in the phonon spectrum. 

Apart from this, there is a possibility that the neigh- 
boring styrene molecules may be electronically affecting 
each other. Although the styrene molecules are not chem- 
ically bonded, there still may be some effect of broad- 
ening. In particular, this effect has been theoretically 
addressed in Refi^ where they conclude that it would 
be higher for unoccupied states than occupied states. 
This conclusion is consistent with Ref4^, where based 
on DFT calculations, it has been reported that for oc- 
cupied states upto 80% of the wave-function is localized 
on the aromatic ring of the styrene molecule. Further- 
more, the unoccupied states are also analyzed theoreti- 
cally in Refs4^i^ reaching the same conclusion. It means 
that hybridization along the chain for occupied state is 
small and hence broadening caused by it is expected to 
be small. Another calculation in Refi^ shows similar be- 
havior with additional information that for the energy 
range correponding to tip voltages less than about -I-2V 
[the voltage range discussed in this paper], the coupling 
between the styrene molecules is small and hence they 
attribute the experimental observation of the slope effect 
to this localized charge on styrene molecules in this volt- 
age range. However, for larger energy range and for tip 
voltages greater than about -I-2V [a case not discussed in 
this paper] , the spectral densities on neighboring styrene 
molecules start overlaping and charge density becomes 
delocalized and hence slope effect vanishes. Since in our 
study, we are calculating transport properties for occu- 
pied states in 0-2V tip voltage range, this effect is ex- 
pected to be small in our voltage range of interest and 
hence is not included in our calculations. 



IV. CONCLUSIONS 

We have studied occupied level conduction through a 
styrene molecule in the vicinity of a negatively charged 
DB on a n++ H:Si(001)-(2xl) surface. We put for- 
ward that the DB not only affect conduction through 
styrene electrostatically within approximately 100 A but 
also electronically up to lOA from the DB by introduc- 
ing a near-midgap state in the LDOS of neighboring Si 
atoms. Dephasing is further expected to play a signifi- 
cant role in these experiments. 



V. APPENDIX 

In this section we show the derivation of Eq. 2. The 
scattering inflow function (S™), outflow function (Sg"*) 
and broadening function (Fg) are defined asi^: 



^in.out 



d{huj) 
27r 



'\nuj)SG"-P{E+huj)S+D'''''''"'{huj)SG 



d{nw) 

27r 



D''"'{hLo)S[G'\E+hLu)+GP{E~huj)]S+D'"'{hLu)S[G'\ 



where ^^'(nw) = (iV + l)Do{hLu) and D^^ihw) = 
N Do{huj) are emission and absorption dephasing func- 
tions respectively and are related by £)ab^jjem _ 
exp^—huj/kBT). N is the equilibrium number of phonons 
given by Bose-Einstein statistics as l/[exp{hLLj/kBT) — 1]. 
For huj < ksT, N + 1 ^ N kBT/hw by keeping first 
term in the Taylor series expansion of exp{fiw/kBT) ss 
1 + hw/ksT, which leads to. 



^s{E) 



d{hLu) 
2tt 



[D*="(;ia;) + D'"'{huj)]SA{E)S 



and 



keT 



Do{hiu) 



Substituting these results in the above Eq. gives Eq. 
d{hui) 2kBDo{huj) 



Ts{E) « T 



27r 



■SA{E)S 



D 



The presence of low lying phonon modes depends on 
many details (type of molecule, surface, bias voltage, 
etc). This implies that not only the energy (hio) of the 
phonon modes would change under different conditions, 
but also the associated dephasing functions Dg. Wc leave 
the detailed study of these phonon modes for the styrene 
chain on Si surface in the presence of a dangling bond 
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for future work and in this study, treat this effect in an 
average manner by using the dephasing strength D as an 
adjusting parameter. 
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